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We present a search for tlic standard model Higgs boson in H WW^*^ decays with e+e^, 
e*/i^ and fJ,'^ fi^ final states in pp collisions at a centcr-of-niass energy of ^/s — 1.96 TeV. The data, 
collected from April 2002 to June 2004 with the D0 detector, correspond to an integrated luminosity 
of 300 - 325 pb~^, depending on the final state. The number of events observed is consistent with 
the expectation from backgrounds. Limits from the combination of all three channels on the Higgs 
production cross section times branching ratio a x BR{H — » WW'*') are presented. 

PACS numbers; 13.85.Rm,14.80.Bn 



In the standard model (SM), the hypothetical Higgs boson is crucial to the understanding of electroweak sym- 
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metry breaking (EWSB) and the mass generation of elec- 
troweak gauge bosons and fermions. Spontaneous EWSB 
predicts the existence of this neutral scalar particle with 
mass Mh, a free parameter in the SM. Direct searches at 
the CERN e+e" collider (lep) yield a lower limit for the 
Higgs boson mass oi Mh> 114.4 GeV Q at the 95% CL. 
Indirect measurements via fits to the electroweak preci- 
sion data give an upper bound of Mh < 280 GeV 2\ at 
the 95% CL. 

In this Letter, we present a search for the Higgs bo- 
son inH ^ WW^*^ -> iiyi'v' {£, £' = e, fj., r) decays with 
e^e^, e^/i^, or fJ-^iJ,^ final states. Tau decays are 
detected in their leptonic decay modes to electrons or 
muons. This is the first search for the Higgs boson at a 
hadron collider in this decay channel, which plays an im- 
portant role in the overall discovery potential of the Higgs 
boson at the Fermilab Tevatron Collider Q . We use data 
collected by the D0 detector between April 2002 and 
June 2004 in pp collisions at y/s = 1.96 TeV of the Fer- 
milab Tevatron Collider. The integrated luminosities are 
325 ±21 pb-\ 318 ±21 pb^^ and 299 ±19 ph'^ for the 
e~^e~, e^fi^ and fj,~^fi~ channels, respectively. The differ- 
ences in the integrated luminosities for various channels 
are primarily due to different trigger conditions. Next-to- 
leading order (NLO) calculations 0, 0| predict the prod- 
uct of the SM Higgs boson production cross section and 
the branching ratio a{pp H) x BR{H WW^*^) of 
11-250 fb for the Higgs masses between 100 and 200 GeV. 
The dominant contribution to the cross section comes 
from gluon-gluon fusion. Extensions of the SM includ- 
ing a fourth fcrmion family |^ predict an enhanced Higgs 
boson production cross section. 

We briefly describe the main components of the D0 
Run II detector Q important to this analysis. The 
central tracking system consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both 
located within a 2.0 T axial magnetic field. The SMT 
strips have a typical pitch of 50-80 /J,m, and the design 
is optimized for tracking and vertexing over the pseudo- 
rapidity range \r]\ < 3, where rj = — In (tan |), where 6 
is the polar angle with respect to the proton beam. The 
system has a six-barrel longitudinal structure, each with 
a set of four silicon layers arranged axially around the 
beam pipe, interspersed with sixteen radial disks. The 
CFT with full coverage for < 1.6 has eight thin coax- 
ial barrels, each supporting two doublets of overlapping 
scintillating fibers, one doublet parallel to the beam axis, 
the other alternating by ±3° relative to the beam axis. 

A liquid- argon/uranium calorimeter surrounds the 
central tracking system and consists of a central calorime- 
ter (CC) covering to \r]\ « 1.1, and two end cap calorime- 
ters (EC) extending coverage to {rjl < 4.2, all housed in 
separate cryostats ^J. Scintillators between the CC and 
EC cryostats provide additional sampling of showers for 
1.1 < |?7| < 1.4. 

The muon system is located outside the calorimeters 



and consists of a layer of tracking detectors and scintil- 
lation trigger counters inside iron toroid magnets which 
provide a 1.8 T magnetic field, followed by two similar 
layers behind each toroid. Tracking in the muon system 
for I77I < 1 relies on 10 cm wide drift tubes jS], while 1 cm 
mini-drift tubes are used for 1 < [yyl < 2 

The H — > W^IF'*) candidates are selected by single or 
di-lepton triggers using a three-level trigger system. The 
first trigger level uses hardware to select electron can- 
didates based on energy deposition in the electromag- 
netic part of the calorimeter and selects muon candidates 
formed by hits in two layers of the muon scintillator sys- 
tem. Digital signal processors at the second trigger level 
form muon track candidate segments defined by hits in 
the muon drift chambers and scintillators. At the third 
level, software algorithms running on a computing farm 
and exploiting the full event information are used to make 
the final selection of events which are recorded for offline 
analysis. 

In the offline analysis, electrons are identified as elec- 
tromagnetic showers in the calorimeter. These showers 
are selected by comparing the longitudinal and transverse 
shower profiles to those expected of the electrons. The 
showers must be isolated, deposit most of their energy 
in the electromagnetic part of the calorimeter, and pass 
a likelihood criterion that includes a spatial track match 
and, in the CC region, an E/p requirement, where E is 
the energy of the calorimeter cluster and p is the momen- 
tum of the track. All electrons are required to be in the 
pseudorapidity range I77I < 3.0. The transverse momen- 
tum measurement of the electrons is based on calorimeter 
cell and track information. 

Muon tracks are reconstructed from hits in the wire 
chambers and scintillators in the muon system and must 
match a track in the central tracker. To select isolated 
muons, the scalar sum of the transverse momentum of all 
tracks other than that of the muon in a cone oi TZ = 0.5 
around the muon track must be less than 4 GeV, where 
TZ = (A(/!))2 + (At;)^ and </> is the azimuthal angle. 
Muon detection is performed over the full coverage of 
the muon system \r\\ < 2.0. Muons from cosmic rays are 
rejected by requiring a timing criterion on the hits in the 
scintillator layers as well as applying restrictions on the 
position of the muon track with respect to the primary 
vertex. 

The decay of the two W bosons into electrons or muons 
results in three different final states e~^e~ + X (ee chan- 
nel), e^/i^+A" {efi channel) and ^^/i~+A {/ifi channel), 
each of which consists of two oppositely charged isolated 
leptons with high transverse momentum and large miss- 
ing transverse energy, fj,, due to the undetected neutri- 
nos. The selection criteria for each channel were chosen 
to minimize the cross section upper limit on Higgs pro- 
duction expected in the absence of signal. To take into 
account the signal kinematic characteristics that change 
with the Higgs boson mass, Mh, some selection cuts 
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are Mh dependent ^3 ■ Higgs boson masses from 
100 GeV to 200 GeV have been studied. 

In all three channels, two leptons originating from the 
same vertex are required to be of opposite charge, and 
must have transverse momenta pr > 15 GeV for the lead- 
ing lepton and pr > 10 GeV for the traihng one (Cut 1). 
Figure shows the good agreement between data and 
Monte Carlo simulation (MC) in distributions of the az- 
imuthal opening angle A(j>u/ between the two leptons for 
the ee (a), the /i/i (c) and the efj, channel (e) after apply- 
ing the lepton transverse momentum cuts. 

In all cases, the background is largely dominated by 
Z/7* production which is further suppressed by requiring 
]^rp > 20 GeV in all three channels (Cut 2). Background 
events are also removed if the Ifj, has a large contribution 
from the mis-measurement of jet energy. The fluctuation 
in the measurement of jet energy in the transverse plane 
can be approximated by AE^°^ ■ sin 9^°^ where Ai?j°' is 
proportional to VE^- The opening angle A<j) (jet, l^j,) 
between this projected energy fluctuation and the miss- 
ing transverse momentum provides a measure of the con- 
tribution of the jet to the missing transverse energy. The 
scaled missing transverse energy defined as 



(A£;j^t • sin0j<=t • cos A(/. (jet, $^)y 



(1) 



is required to be greater than 15 (Cut 3). 

The charged lepton system and the neutrinos are emit- 
ted mostly back-to-back, so the invariant mass for the 
leptons from the Higgs decay is restricted to Majl. 
Thus, the invariant mass is required to be mu < 
Mh/2 (Cut 4). In the ee channel the cut is altered to 
niee < TOm(80GeV, Mh/2). In the /i/i channel a lower 
cut boundary with m^^ > 20 GeV is required to remove 
events from J/ip, T and Z/j* production. The sum of 
the pt of the leptons and l^j, is required to be in the range 
Mh/2 + 20 GeV < p^! + p^^ + < Mh for the ee and 
e/i channel and Mh /2 + 10 GeV < + + fj, < Mh 
for the fifi channel (Cut 5). The transverse mass, defined 

as to"' = y^2p"'^j,(l — cos A(/)(p"', with the di- 

lepton transverse momentum pf^ , should be in the range 
Mh/2 < m^rf' < Mh - 10 GeV (Cut 6). The latter two 
cuts reject events from W^-f jet/7 ^-nd WW production 
and further reduce backgrounds from Z/j* production. 
Finally, to suppress the background from tt production, 
the scalar sum of the transverse energies of all jets with 
> 20 GeV and \r]\ < 2.5, Ht, is required to be less 
than 100 GeV (Cut 7). Remaining Z boson and multi-jet 
events can be rejected with a cut on the opening angle, 
A<j)u' < 2.0 (Cut 8), since most of the backgrounds ex- 
hibit a back-to-back topology. This is not the case for 
Higgs boson decays because of the spin correlations in the 
decay. Figure shows the distributions of the azimuthal 
opening angle A(j)u' between the two leptons for the ee 
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FIG. 1: Distribution of the opening angle Acpui after ap- 
plying the initial transverse momentum cuts in the (a) ee, 
(c) fifj, and (e) e/i channel. Figures (b), (d) and (f) show 
the A(j>ee' distributions after the final selection except for the 
A<j>iii criterion for the ee, /t/x, and e/i channel, respectively. 
The arrows indicate the cut values. The QCD contribution is 
negligible in Figs, (c) and (d). 



(b), the fj,fi (d) and the e/z channel (f) before applying 
the final cut on Acjjui . 

To maximize the sensitivity, the selection in the fifi 
channel is slightly changed for Higgs boson masses Mh= 
140 and 160 GeV. For a better Z/j* background sup- 
pression cuts 4, 5 and 6 are replaced by the following 
cuts: the invariant mass to^^ should be in the range 
20 GeV < m^^ < 80 GeV (Cut 4). Since the momentum 
resolution is degraded for high pr tracks, an additional 
constrained fit is performed to reject events compatible 
with Z boson production (Cut 5). The sum of the muon 
transverse momenta and the missing transverse energy 
should be + pl^ + > 90 GeV (Cut 6). 

The efficiency for H —>■ WW^*'' iv £'v' signal events 
to pass the acceptance and selection criteria is deter- 
mined using the pythia 6.2 111! event generator followed 
by a detailed GEANT-based 12] simulation of the D0 de- 
tector. All trigger, reconstruction and identification effi- 
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TABLE I: Overall detection efficiencies (in %) for 
H WW^*^ Iv i'u' decays for the three channels after all 
cuts. Quoted are the overall uncertainties, combining statis- 
tical and systematic components in quadrature. 

Mh (GeV) ee e/i /i/i 

Too 0.56 ±0.05 1.02 ± 0.06 0.44 ± 0.03 

120 1.18 ±0.09 2.0 ±0.1 1.02 ±0.06 

140 1.55 ±0.08 2.9 ±0.2 1.34 ±0.08 

160 2.1 ±0.1 3.9 ±0.2 2.0 ±0.1 

180 2.1 ±0.1 3.9 ±0.2 1.68 ±0.09 

200 1.57 ±0.09 3.2 ±0.1 1.53 ±0.07 



ciencies are derived from the data. The kinematic accep- 
tance efficiency is derived from MC. The overall detec- 
tion efRciencies range from (0.44±0.03)% to (3.9±0.2)% 
depending on the decay channel and Mh- Tabic sum- 
marizes these efficiencies. 

Using the NLO cross sections calculated with higlu 
and HDECAY and the branching ratio BR of 0.1068 ± 
0.0012 for W ^ iiy the expected number of events 
for H VFW^^*^ decays from all three channels is 0.68 ± 
±0.03 (syst) ± 0.04 (lum) events for a Higgs boson mass 
Mh= 160 GeV. The signal expectation for different Higgs 
masses Mh are given in Table ITU 

Background contributions from Z/^f*, VF±jet/7, ti, 
WW, WZ and ZZ events are estimated using the 
PYTHIA event generator normalized to their NLO cross 
sections In addition, W+jet/^ contributions are ver- 
ified using ALPGEN Ti'l . All events are processed through 
the full detector simulation. The background due to 
multi-jet production, when a jet is misidentified as an 
electron, is determined from the data using a sample 
of like-sign di-lepton events with inverted lepton quality 
cuts (called "QCD" in Fig. A summary of the back- 
ground contributions together with signal expectations 
and events observed in the data after the final selection 
is shown in Table ITll There is good agreement between 
the number of events observed in the data and the various 
backgrounds in all three channels. The largest difference 
between the data and the background expectation, at 
Mh= 120 GeV, corresponds to a background probabil- 
ity of 6%. The e/i channel has both the highest signal 
efficiency and best signal-to-background ratio. 

Various sources of systematic uncertainties that affect 
the background estimation and signal efficiencies have 
been studied. In these calculations, parameters are var- 
ied within ±1(7 of their nominal values, where a is deter- 
mined by the corresponding uncertainties. The trigger 
efficiency, electron and muon identification efficiencies, 
jet energy scale, electron and muon momentum resolu- 
tion, parton distribution function uncertainty and cross 
section calculation of Z/ 7* , WW and tt events contribute 
to the systematic uncertainties. The total systematic un- 
certainties for the background estimate and signal effi- 
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FIG. 2: Expected and observed upper limits on the cross sec- 
tion times branching ratio a x BR{H ~* WW'*') at the 95% 
CL together with expectations from standard model Higgs bo- 
son production and an alternative model. The lep limit on 
the standard model Higgs boson production is taken from Q 
and the 4th generation model prediction is described in 1^. 

ciencies for the six Higgs boson masses are given in Ta- 
blc lIIII The largest contribution to the systematic uncer- 
tainty on the background for small Higgs boson masses 
comes from the jet energy scale due to the large M^±jet/7 
background, whereas for high Higgs boson masses the W 
boson pair production cross section gives the largest sys- 
tematic uncertainty. The uncertainty of the parton dis- 
tribution function is the largest uncertainty of the signal 
efficiency. The uncertainty of the luminosity measure- 
ment is 6.5%. 

Since the remaining candidate events after the selec- 
tion are consistent with the background expectation, lim- 
its on the production cross section times branching ratio 
(T X BR{H WW^*^) are derived using a modified fre- 
quentist method described in Ref. . It provides the 
confidence level for the background to represent the data, 
CLb, and the confidence level for the sum of signal and 
background hypothesis CLs+b- The 95% CL limit is 
obtained by requiring CLs+b/GLb — 0.05. The uncer- 
tainty on the background and the expected signal events 
were determined from the statistical and systematic un- 
certainties and luminosity uncertainty. Table IIVI shows 
the individual expected and observed upper limits on the 
cross section times branching ratio a x 
for the combination of the three different decay channels 
for six different Higgs boson masses. The different val- 
ues of the upper limits are due to different background 
expectations and signal efficiencies. The best limits are 
achieved for large Higgs masses since background expec- 
tations decrease while signal efficiencies increase. 

FigureElshows the expected and observed cross section 
limits for a x BR{H WW^*^) for the different Higgs 
boson masses compared with predictions from the SM 
and from an extension including a fourth fermion fam- 
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TABLE II; Number of signal and background events expected and number of events observed after all selections are applied. 
Only statistical uncertainties are given. 



Mh (GeV) 


100 


120 


140 


160 


180 


200 




0.007 ± 0.001 


0.125 ± 0.002 


0.398 ± 0.008 


0.68 ± 0.01 


0.463 ± 0.009 


0.210 ± 0.004 


zh' 


7.9 ± 1.1 


7.5 ± 1.0 


3.8 ±0.6 


4.0 ±0.7 


6.6 ±0.9 


9.9 ± 1.1 


Diboson 


4.4 ±0.2 


8.1 ± 0.2 


11.7 ±0.3 


12.3 ±0.3 


11.6 ±0.3 


9.6 ±0.3 


ti 


0.03 ± 0.01 


0.11 ± 0.02 


0.29 ±0.02 


0.47 ±0.03 


0.66 ± 0.05 


0.72 ± 0.05 




16.9 ±2.2 


14.2 ±2.1 


5.8 ± 1.2 


2.8 ±0.9 


0.7 ±0.5 


0.7 ±0.5 


Multi-jet 


0.6 ±0.3 


0.3 ±0.1 


0.2 ±0.1 


0.2 ±0.1 


0.3 ±0.1 


0.3 ±0.1 


Background sum 


29.9 ±2.5 


30.1 ± 2.3 


21.8 ± 1.4 


19.7 ± 1.2 


19.8 ± 1.1 


21.2 ± 1.2 


Data 


27 


21 


20 


19 


19 


14 



TABLE 111: Systematic uncertainties (in %) of the 
H WW'^*^ signal efficiencies and of the number of back- 
ground events (BG), for the ee, ejj, and fijj, channels. Uncer- 
tainties exclude the uncertainty of the luminosity measure- 
ment. 



ee e/i /i/i 



M„{GeV) 


Signal 


BG 


Signal 


BG Signal BG 


100 


8.3 


9.5 


6.4 


11.4 7.8 7.2 


120 


8.3 


8.6 


6.7 


13.6 7.3 7.5 


140 


6.4 


6.7 


6.9 


13.6 7.2 8.3 


160 


6.6 


7.3 


6.7 


12.0 7.1 8.3 


180 


6.9 


10.3 


6.6 


13.0 7.3 14.6 


200 


6.8 


10.6 


6.1 


12.3 6.9 18.1 


TABLE IV: 


Expected and observed upper limits at the 95% 


CL for the Higgs boson production cross section times branch- 


ing ratio a X 


BR{H ^ 


WW 


*') for six values of Mh- 


Mh (GeV) 




100 


120 


140 160 180 200 


Expected limits (pb) 


20.3 


9.5 


5.9 4.0 3.9 4.5 


Observed limits (pb) 


18.5 


5.6 


4.9 3.7 4.1 3.2 



ily 1^1 . With the current dataset, no region of the SM 
prediction can be excluded. 

To conclude, we have searched for the Higgs boson in 
H ^WW^*'> ^Ivt'v' {£,£' ^e,n,T) decays with e+e", 
e^^^ or ^+^~final states in pp collisions at y/s — 1.96 
TeV. The data is consistent with the expectation from 
backgrounds. Since no excess has been observed, limits 
on the production cross section times branching ratio a x 
BR{H VFVF(*)) have been derived. 
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